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Generally, the conditions for deep sub-Doppler laser cooling do not match the conditions for 
the strong atomic localization that takes a place in deeper optical potential and, in consequence, 
leads to larger temperature. Moreover, for a given detuning in a deep optical potential the secular 
approximation which is usually used for quantum description of laser cooling becomes no more valid. 
Here we perform an analysis of atomic localization in optical potential based on a full quantum 
approach for atomic density matrix. We also show that the laser cooling in a deep far-off detuned 
optical potential, created by a light field with a polarization gradient, can be used as an alternative 
method for forming high contrast spatially localized structures of atoms for the purposes of atom 
lithography and atomic nanofabrication. Finally, we perform an analysis of the possible limits for 
the width and the contrast of localized atomic structures that can in principle be reached by this 
type of the light mask. 

PACS numbers: 32.80.Pj, 42.50.Vk 



I. INTRODUCTION 

Laser cooling and manipulation of neutral atoms is one 
of the priority field of atom optics. In the recent time the 
major developments and success have been obtained in 
atom lithography and direct deposition of atoms utiliz- 
ing lirfit fields as an immaterial optical mask for atomic 
beam In most nanofabrication experiments, atomic 

structures are realized by far off detuned periodical con- 
servative potential created by intense laser fields acting 
as an array of immaterial light lenses for the atoms. The 
influence of the spontaneous emission on the focusing is 
considered to be negligible because of the large light de- 
tuning and short interaction times. In essence, the atom 
trajectory affected by conservative dipole force without 
any loses (or dissipation) of energy in the atomic beam. 
In this case the atomic beam focusing has a classical anal- 
ogy and can be described with methods developed for par- 
ticle optics Q. As in normal optics, the feature size is 
limited by a combination of chromatic aberration caused 
by the broad longitudinal and transverse velocity distri- 
bution of an atomic beam. Therefore an additional laser 
cooling field is required to prepare a well collimated and 
transversely cooled atomic beam to minimize deleteri- 
ous effects. Additionally, because of spherical aberration 
some atoms do not focus well and contribute to pedestal 
background. These factors are dominant and do not al- 
low one to reach the theoretically predicted diffraction 
limit for atom optics determined by de Broglie wave- 
length of atoms, only a few picometres. Therefore the 
new alternative methods for atom lithography are inten- 
sively investigated. 

Recently, the idea of combining the traditional focus- 
ing method with the well known concept of laser cool- 
ing was suggested for a blue detuned intensive light field 



H, 0] and was mentioned earlier in Q. Here the in- 
tense light field, commonly used for focusing, was used 
to create a deep optical potential. Additional dissipa- 
tive light force provides additional cooling of the atoms 
to the minimum of optical potential at blue detuning. 
The characteristic time when the dissipation processes 
take effect is a few inverse recoil frequencies uj^ 1 (where 
hujji = h 2 k 2 /2M is recoil energy, gained by an atom with 
mass M at rest after emission of a photon with momen- 
tum hk). This time is several tens of microseconds for 
the number of elements with closed dipole optical transi- 
tions suitable for laser cooling. Thus for commonly used 
atomic beams with thermal longitudinal velocity it might 
be difficult to realize this type of dissiaptive optical mask 
experimentally due to power limitations of laser system 
used. 

In the following we consider the alternative regime of 
dissipative optical mask, created by red detuned low in- 
tensity light field with nonuniform polarization. It is 
well known that low intensity light field with polariza- 
tion gradients can be used for sub-Doppler laser cooling 
of neutral atoms. This mechanism of laser cooling is well 
understood and has been thoroughly studied by a num- 
ber of authors 0, @, 0, IT^. especially with respect 
to the temperature of laser cooling, the atomic momen- 
tum distribution, and the localization IWL llfll that can 
be measured by spectroscopy methods [171 llq. Due to 
the extremely complex master equation for the quantum 
description of atomic motion in light fields the secular 
approximation 0, ITllfT^ . \T?l \Tty was initially suggested 
is valid in the limit [10] 

ViVfi^« \s\h- (1) 

This limit assumes the energy separation between differ- 
ent energy bands of atoms in optical potential are much 
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greater then their width due to optical pumping and tun- 
neling effects. Here the light-shift well depth Uq define 
the optical potential depth. 6 = uj — ui is the detuning 
between the laser uj and atomic transition ujq frequencies, 
and 7 is the radiative decay rate. This approximation is 
valid, for a given potential depth, in the limit of large 
detuning. On the contrary, it might failed in a deep po- 
tential at given detuning. More over, even if the secular 
approximation is well fulfilled for the lowest vibrational 
levels it might break down for upper ones where the sep- 
aration between energy bands becomes smaller due to 
potential anharmonicity effects, and especially for above- 
barrier motioned atoms. 

In the present work we investigate the applicability 
of laser cooling in a deep optical potential created by 
light field with nonuniform polarization for generating 
spatially localized atom structures with high contrast for 
atom lithography. We consider the conditions far from 
the situation of extremely low sub-Doppler cooling cases. 
Thus, to describe the localization of atoms more correctly 
in here we do not restrict our consideration by secu- 
lar approximation. Rather we perform a full quantum 
numerical analysis of generalized optical Bloch equation 
for atomic density matrix elements. In particular, we 
consider the light field parameters beyond the secular 
approximation limit. Finally we analyze the structure 
width and contrast of localized atoms, important param- 
eters for technological applications. 



II. MASTER EQUATIONS 

Let us consider one-dimensional (along the z axis) mo- 
tion of atoms with total angular momenta j g in the 
ground state and j e in the excited state in a field of two 
oppositely propagating waves with the same frequency 
and intensity 



E(z,t) = E (ei e ikz + e 2 e 
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elements [llll^: 



dt P 



f{p} 



where the Hamiltonian H is given by 
H = ^ + H5SV*V. 



(4) 
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The last term in (0 describes atoms interaction with the 
light field in resonance approximation, where 



V = Vi e + V2 e 
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and operator T q is written through the Clebsch-Gordan 
coefficients: 



Ms I 
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written in the basis of sublevel wave functions for exited 
\j e ,fi e ) and the ground \j g ,n g ) atomic states. 

In addition, the relaxation part of kinetic equation for 
atomic density matrix <@J has the following form 



f{p} = -^{vty,p} 



(8) 
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ft e -iksi v e -^f q K q {s) ds 



where {a, c} = dc + cd standard anticommutator defi- 
nition and z is position operator. This term describes 
redistribution of the atom on the ground state energy 
sublevels with taking into account the recoil effects in 
spontaneous photon emission. The functions K± 1 (s) = 
(1 + s 2 )3/8 and K (s) = (1 - s 2 )3/4 is determined by 
the probability of emission of a photon with polarization 
q = ±1,0 into direction s = cos(6>) (relative to the z 
axis) . 



Here E is the amplitude of each of the oppositely prop- 
agating waves. The unit vectors ei and e 2 determine 
their polarizations with components in cyclic basis 
{e = e z ,e ±1 = =F(e^ ±ie y )/V2}. 

In this work we restrict our consideretion by the weak- 
field limit, i.e. with small saturation parameter 
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n 2 



(5 2 + 7 2 /4 



(3) 



Here fl = —Egd/H is the single-beam Rabi frequency 
characterize the coupling between the atomic dipole d 
and light field. 

In the weak-field limit, the atomic exited-state can be 
adiabatically eliminated, and atom motion is describing 
by reduced equation for the ground-state density matrix 



III. EQUILIBRIUM ATOMIC DENSITY 
MATRIX 



There are number of approaches developed for calcu- 
lation of the evolution of the atomic density matrix. The 
quantum problem is more difficult because it incorpo- 
rate evolution of numbers of internal and external com- 
ponents of the density matrix. The majority of works 
are based on secular approximation for density matrix 
elements [T3. fiH. H^|. It consists in the following: First, 
the eigenstates and eigenenergies of the Hamiltonian H 
should be found. Then we can consider only the evolu- 
tion of diagonal elements of atomic density matrix in this 
eigenstatas basis. It can be written in a form of balance 
equation with the rates characterizing the relaxation part 
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of the master equation. The secular approximation for 
the general master equation <@J is valid for the lowest 
vibrational energy levels when the the energy separation 
between different bands are much greater than the ef- 
fective width due to optical width and tunneling effects. 
This condition implies very large detuning. However, the 
density of energy states increases for the upper energy 
levels and the energy difference between adjacent states 
could be very small This circumstances therefore 

make very hard to use the secular approximation for de- 
scribing of hot and nonlocalized atomic fraction. 

In order to take into account effects of localization in 
optical potential and modulation depth of the spatial dis- 
tribution of atoms more correctly we utilize a different 
approach for the quantum calculation of the master equa- 
tion JUl. 

In the Wigner representation for atomic density matrix 
p(x,p) the general master equation <@J takes the follow- 
ing form: 

lp = -i6s[v^V,p]-^-{v^V,p}+m (9) 

where commutator of density matrix with kinetic part of 
Hamiltonian (jj) reduces to partial deviation on z 

The field part of Hamiltonian © has only the zeroth and 
the second spatial harmonics: 



V*V = Wo + W+e l2kz + W-e 



i2kz 



(11) 



thus the commutator and anticommutator in the right 
hand side of equation (jHJ could be written as: 



V'V p =F P V^V = W OP (z,p)Tp(z,P)Wo 

+ (\Y- p(z,p + hk)=f p(z,p- hk)W- 

+ (w + p(z,p-hk)TP{z,P + tik)W + 



(12) 

p — i2kz 
Alhz 



The last term describing relaxation due to spontaneous 
emission of photons has well-known form in Wigner rep- 
resentation: 

/hk 
dp'/hk K q ( P '/hk) ft 

Vi p{z,p + p') V 2 f e l2kz + V 2 p{z, P + p') V?e- i2kz 



+ v lp (z,p + p' - hk) + v 2 p(z, p + P ' + hk) y 2 1 



Equation © admits solution that is periodic in the posi- 
tion variable. We further introduce a Fourier-series ex- 
pantion for atomic density matrix in spatial coordinates 



«*,p)=5> (n) (p) 
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FIG. 1: Steady-state spatial (a) and momentum (b) distribu- 
tions, and the total population of spatial harmonics of atomic 
ground-state density matrix (c) for atoms with j g = 1 — » j e = 
2 optical transition and mass as for Cr atoms. The light field 
detuning 8 = — 4O7 and saturation parameter S = 0.5. 



and rewrite the master equation for discrete Fourier com- 
ponents of density matrix p^ : 



d_ 



AI 



(14) 



{/><«> }+C + {p^} + C^{p^} . 



In steady-state problem (d/dt p = 0) such a recursion 
may often be solved by the method of continued fraction. 
This approach is used for solution of the optical-Bloch 
equations in different spectroscopy tasks as well as for 
the calculation of the force on atom in the light field (see 
for example 0,H3]). The major distinction here is that 
these equations for density matrix contain the recoil ef- 
fects that makes them more complicated for calculation. 
Additionally we note, the similar approach was described 
in |2ll where the authors analyzed the laser cooling (ve- 
locity distribution) of two-level atom in the recoil limit 
and thus restrict their consideration only by the zeroth 
spatial harmonics for the ground state density matrix. 
In our case the number of the spatial harmonics depends 
on the light field and atomic parameters. Typically we 
use less than 30 harmonics that is enough to obtain the 
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Element 


cooling 


M 


A 


Is 




transition 




(nm) 


(mW/cm 2 ) 


7 Li 


2 2 Si/2 — * 2 2 P 3 / 2 


46 


671 


2.56 


23 Na 


3 01/2 — > 3 P3/2 


198 


589 


6.34 


39 R 


4 2 Si/ 2 -» 4 2 P 3/2 


358 


766 


1.81 


85 Rb 


5 S1/2 — > 5 P3/2 


770 


780 


1.63 


133 Cs 


6 2 5 1/2 -> 6 2 P 3/2 


1270 


852.3 


1.06 


52 Cr 


4 7 S , 3 -> 4 7 P 4 


115 


425.6 


8.49 


27 Al 


3p 2 P 3/2 3d 2 P 5/2 


85 


309.4 


57 


69 Ga 


4p 2 P 3 /2 -» 4d 2 £> 6/2 


382 


294.4 


127 


115 In 


5p 2 P 3/2 -> 5d 2 P 5/2 


634 


325.7 


78 


107 Ag 


5 S1/2 — * 5 P3/2 


601 


328 


76.8 



TABLE I: Dimmensionless atomic mass parameter corre- 
sponding to laser cooling transition for different elements. 



spatial solution for equilibrium atomic density matrix in 
considered range of light field parameters. 

The spatial and momentum steady-state distribution 
of atoms with j g = 1 — > j e = 2 optical transition for 
S = — 4O7, S = 0.5 and chromium mass are shown 
in FigQIa) and (b). The Fig^c) represents the to- 
tal population of spatial harmonics of atomic ground- 
state density matrix integrated over momentum space 
= J Tr{p^ n '(p)} dp. The zeroth harmonic is equal 
to 1 that is the normalization condition. As it seen here 
the population of higher harmonics are rapidly decrease 
with number n. 



IV. RESULTS 

In this section we turn our attention to the steady-state 
spatial distribution of the atoms in the optical potential 
created by the light field with lin _L lin configuration. 
We choose this configuration as a brightest example of 
the light field with nonuniform polarization. Only light 
field ellipticity varies in position space while the other 
parameters (intensity, phase, orientation of polarization 
vector) stay unchanged. More over, the optical potential 
created by this field configuration has a period of A/4 
that makes it very attractive for deposition of atomic 
structures with high spatial periodicity. 

There are several physical parameters that are required 
to characterize a given laser cooling situation. These are 
the atomic mass M, the wavelength A, and the natu- 
ral line-width 7. In addition, there are two light-field 
parameters in the low-field limit: detuning S, and sat- 
uration parameter S. It is possible to choose reduced 
dimensionless units (ft = 1, k = 1, 7 = 1), so that the 
dimensionless atomic mass M can be defined from rela- 
tion 7/wr = 2M @. This is so called a quasiclassical 
parameter that characterizes an atom's kinetic in a light 
field. In particular, it describes the rate of kinetic pro- 
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FIG. 2: Spatial FWHM width of localized atomic structures 
(a), and contrast (b) as function of parameter 8 /i for different 
field detuning (S/-j = -5, -10, -20, -40, -80, -160 from top 
to bottom in (a) and from bottom to top in (b)) in a model of 
atom with j g = 1/2 —* j e =3/2 optical transition. 



cesses and evolution of atomic distribution function in 
momentum space, thus the typical cooling time is order 
of r = (jJr/(jS). 

We perform an analysis of the effects of atom local- 
ization in optical potential. Contrast can be defined as 
the ratio of spatial distribution modulation depth to its 
amplitude C = h/H (FigQ] (a)) as in 1]. As it can be 
seen directly from ijjjj , there are only two parameters that 
characterize the stationary solution for an atomic density 
matrix. This is the detuning, that can be measured in 
the units of natural line- width S — S/j and the second 
one is non-dimensional parameter /u: 



H = SM . 



(15) 



This scale parameter makes it possible to implement the 
results of calculation for the elements with allowed closed 
dipole optical transition having degeneration over the an- 
gular momentum energy levels (table I). 

Note, that in the secular approximation the sta- 
tionary solution is characterized only by the ratio of op- 
tical potential depth to recoil energy Uq/ojr that is pro- 
portional to 5[i in our notations. Thus we first express 
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FIG. 3: Asymptotic FWHM width of localized atomic struc- 
tures as function of light field detuning in a model of atom 
with j g = 1/2 — » j e = 3/2 optical transition. 



the results in parameters 5 and <5/i. As it seen in the 
Fig. [21 the differences between the curves become more 
significant with the optical potential depth Uo increas- 
ing, corresponding to out of framework from secular ap- 
proach. The width of the localized structures and the 
contrast tends to an asymptotic constant value with in- 
creasing of the light field intensity, however that depends 
on light field detuning. In order to find the asymptotic 
values for the width we fit curves with an empirical low 
w(x) = a + b/^/1 + ex. The results for asymptotic values 
a are shown in Fig. 



The Fig. 0] represents the FWHM width and con- 
trast of the spatial structures as a function of parame- 
ter [i at different detunings for the atomic model with 
j g = 1 — * j e = 2 optical transition. The FWHM width of 
the spatial structures is monotone decrease for a large de- 
tuning and parameter \i. In spite of the fact that the equi- 
librium temperature is growing with the depth of optical 
potential, the localization of the atoms becomes stronger 
with growth of these parameters. Additionally the con- 
trast tends towards its maximum value Fig. EJb). The 
dashed vertical lines on the Fig^l underline the limita- 
tions of the weak-field theory for different elements from 
the table I in assumption S < 0.5. Note, that this is qual- 
itative assumption. For thorough analysis of the weak- 
field theory limits the solution of the quantum equation 
for the total atomic density matrix is required with tak- 
ing into account exited state population. However, the 
width and contrast curves FigQJhave a very strong de- 
pendence on fi, thus the localization effects remain rather 
significant for enough small saturation parameter espe- 
cially for "heavy" atoms table I. 



V. CONCLUSION 

We performed the fully quantum analysis of atomic 
localization obtained by laser cooling in nonuniformly 
polarized low-intensity light field. Generally, conditions 
for a deep laser cooling mismatch conditions for strong 
atomic localization that require deeper optical potential 
and consequently leads to higher laser cooling tempera- 
ture. Additionally, in a deep optical potential the secular 
approximation restricted by relation on light field de- 
tuning and the depth of optical potential. In our treat- 
ment we had no such limitation allowing a more cor- 
rectly describe the spatial solution for atomic distribution 
function including localized and above-barrier motioned 
atoms. The stationary solution of atoms is a function 
of the light field detuning 5 and the dimensionless pa- 
rameter [i dSJ. We analyzed the width and contrast 
of localized atomic structures as a function of these pa- 
rameters. We showed that the atomic structure width 
and contrast have a strong dependence on n and tend to 
constant values with an increasing optical potential dept 
that depends on the light field detuning. 

We demonstrated the applicability of laser cooling in 
far-off detuned deep optical potential, created by a light 
field with polarization gradient, as a dissipative optical 
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FIG. 4: Spatial FWHM width of localized atomic structures 
(a), and contrast (b) as function of parameter fi for different 
fi,eld detuning in a model of atom with j g = 1 — > j e = 2 optical 
transition. 
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mask for the purposes of atom lithography and nanofab- 
rication for generation of spatially localized atomic fea- 
tures with high contrast. This type of the light mask 
can be an alternative method for creation of spatially 
localized atomic structures. The remarkable distinction 
of this method from previous non-dissipative light mask 
is that the suggested one is not sensitive to any aber- 
ration effects. More over, this type of optical mask has 
no classical analog and can not be described by meth- 
ods for classical optics. The width and the contrast of 
localized atomic structures here is determined by atomic 
energy dissipation mechanisms in the light field. Finally, 



we perform analysis of the possible limits for the width 
and the contrast of localized atomic structures that can 
in principle be reached by this type of the light mask. 
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